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Rotationally resolved infrared spectra are reported for the X-HCN (X = CI, Br, I) binary com- 
plexes solvated in helium nanodroplets. These results are directly compared with that obtained 
previously for the corresponding X-HF complexes [J. M. Merritt, J. Kiipper, and R. E. Miller, 
PCCP, 7, 67 (2005)]. For bromine and iodine atoms complexed with HCN, two linear structures 
are observed and assigned to the 2 E 1/ / 2 and 2 n 3 / 2 ground electronic states of the nitrogen and hy- 
drogen bound geometries, respectively. Experiments for HCN + chlorine atoms give rise to only 
a single band which is attributed to the nitrogen bound isomer. That the hydrogen bound isomer 
is not stabilized is rationalized in terms of a lowering of the isomerization barrier by spin-orbit 
coupling. Theoretical calculations with and without spin-orbit coupling have also been performed 
and are compared with our experimental results. The possibility of stabilizing high-energy struc- 
tures containing multiple radicals is discussed, motivated by preliminary spectroscopic evidence for 
the di-radical Br-HCCCN-Br complex. Spectra for the corresponding molecular halogen HCN-X2 
complexes are also presented. 



INTRODUCTION 

Much of our knowledge of intermolecular forces has 
come from the high-resolution spectroscopy of weakly 
bound closed shell complexes Q, and there is consid- 
erable promise that similar advances can be made for 
open shell complexes. Experimentally, the study of open 
shell reactive molecular complexes is more difficult than 
their closed shell counterparts due to the fact that their 
tendency to react must also be surpressed. The cooling 
provided by free jet expansion |2j, La, Li an d matrix isola- 
tion [a,Sl3 is sometimes sufficient to stabilize complexes 
of this type. Recently, nano-scale liquid helium droplets 
have also emerged as a nearly ideal spectroscopic matrix 
for the study of highly metastable species H H Lul El , 
including the stabilization of pre- and post-reactive clus- 
ter systems 0, El Q] . Such weakly bound complexes 
sample the long range van der Waals forces in the en- 
trance and exit channel regions of the potential energy 
surface, which have recently been shown to strongly in- 
fluence the reaction dynamics of CI + HD 0, Q, E} 
and the near threshold photo-dissociation of formalde- 
hyde 0, El- Despite the fact that the magnitude of 
such orientational forces are usually very small in com- 
parison with the energy of the transition state, the torque 
on the reactants at long range acts to deflect trajectories 
towards or away from the transition state [lfij ] . 

"Heavy+light-heavy" systems such as X-HY, where X 
and Y are both halogen atoms, have emerged as proto- 
typical systems to study chemical reaction dynamics at a 
fully quantum mechanical level [13, El El El ■ The 
attraction of such complexes is due to the fact that they 
contain only 2 heavy atoms, thus simplifying the nuclear 
degrees of freedom, yet the corresponding electronic de- 
grees of freedom can be quite complex. Indeed, because 
the unpaired electron of the free halogen atom resides in 



a p orbital, three potential energy surfaces are needed to 
fully characterize the interactions, which correspond to 
the three relative orientations of the p orbital with re- 
spect to HX. These surfaces become degenerate for large 
separations, and the coupling between the surfaces can- 
not be neglected EH- Furthermore, spin-orbit cou pling: 
will also reshape the potential energy surfaces |2f)ll2tll27j |. 
In this work we extend these studies by looking at a re- 
lated X-HY complex, where Y is now the pseudo-halogen 
CN. 

The reactions of halogen atoms with HCN have been 
the focus of extensive exp erimental [28, 29, 3Qj and the- 
oretical study |3, El El| j aimed at elucidating the role 
of the potential energy surface in determining the associ- 
ated reaction dynamics. The CN group is often regarded 
as a pseudo-halogen, and thus there is interest in com- 
paring the reaction dynamics of this tetra-atomic sys- 
tem with the analogous X-HY (X & Y = halogen atoms) 
triatomic systems. Interestingly, the long range interac- 
tions between XH + CN are expected to differ signifi- 
cantly from X + HCN due to the threefold degeneracy 
associated with the 2 P ground state of halogen atoms, 
compared to the nondegenerate 2 S ground state of CN 
radicals, further lowering the symmetry of this hydro- 
gen exchange reaction. Despite this, HCN is an ideal 
test case because the CN group is usually found to act 
as a single unit, but there is still considerable debate 
about whether or not the CN bond is a spectator to the 
dynamics of many different reactions El E3. EH ■ The 
additional degrees of freedom that the HCN molecule 
presents, open up many new channels of the potential 
energy surface that can be explored, and the X + HCN 
systems have provided benchmarks for the theoretical 
treatment of such multi-dimensional dynamics. 

The Cl-H-C-N potential energy surface has been stud- 
ied experimentally from three different starting positions. 
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De Juan et al. 33] have investigated the H + C1CN — > 
HC1 + CN reaction by colliding translationally hot H 
atoms formed from the 248 nm photolysis of H 2 S, with 
C1CN. This experimental arrangement prepares the re- 
actants with 21.6 kcal mol -1 of collision energy, which is 
~43 kcal mol^ 1 above the energy of CI + HCN (see Fig- 
ure^. By analyzing the nascent CN rovibrational popu- 
lations, as well as the Doppler profiles of the vibrational 
transitions, it was proposed that CN is a spectator under 
these conditions, and that most of the available energy 
goes into HC1 vibration and into translational energy of 
the two fragments. Sims and Smith [34| and Frost et al. 
[13 have studied the RG\(v HC i) + CN(v C n) -> HCN + 
CI reaction, and came to the conclusion that CN is also a 
spectator in this reaction, due to the fact that there is a 
negligible enhancement of the reaction rate upon exciting 
the CN stretching vibration. A third set of experiments 
were carried out by Metz et al. 

mm 

and Krcher et 

al. focusing on the reaction of translationally hot 

chlorine atoms with highly vibrationally excited HCN: 

Cl + HCN{v u v 2 ,v 3 ) -> HCl + CN 

Metz and coworkers showed that HCN, when prepared 
with either 4 quanta of CH stretch (004) or 3 quanta of 
CN stretch plus two quanta of CH stretch (302) react 
with comparable rates, and conclude that under these 
conditions CN "is clearly not a spectator" [8j| . An- 
other interesting observation from this work is that only 
14% of the available energy was found as HC1 vibration, 
which is surprisingly small since this is the bond that is 
formed in the reaction. To justify this observation an 
addition-elimination mechanism was proposed which in- 
volved an intermediate HC1CN product, as opposed to 
direct collinear abstraction. Indeed, experiments in ar- 
gon matrices have now proven the existence of HFCN 
|38l m| . If the intermediate survived long enough in the 
gas-phase to permit energy redistribution, this could ac- 
count for the transfer of energy from HCl(u) to CN(i>). 
According to this mechanism the CN cannot be a simple 
spectator. Kreher et al. confirm these observations, how- 
ever point out that HCN(302) does give rise to higher vi- 
brational excitation of CN than does HCN (004), so some 
memory of the initial state is retained. 

The global potential energy surface for CI + HCN 
was first studied theoretically by de Juan et al. at the 
MP4/3-21G* level to aid in the interpretation of their 
experimental results j^- Their calculations did show 
the existence of a covalently bound HC1CN intermediate 
product, but the relative energies of the reactants and 
products showed significant discrepancies with the exper- 
imentally determined values. Harding later revisited the 
Cl-HCN surface at the RHF+l+2+QC/cc-pVDZ level 
and found the energies to be in much better agreement 
with experiment [22|. Although the overall topology of 
the two surfaces are similar, the existence of a direct 
HC1CN — > HC1 + CN pathway originally proposed by de 



Juan et al. was not found by Harding. Instead, Harding 
points out that they most likely found the transition state 
for collinear abstraction. Currently an exhaustive search 
of the six-dimensional potential energy surface is not pos- 
sible, however Harding states that this HC1 elimination 
pathway is unlikely based on the overall topology of the 
potential. He finds that the transition state of the CI 
+ HCN — » HC1 + CN reaction is a collinear abstraction 
mechanism, and that the reaction path does not directly 
involve the HC1CN complex, in contrast to that proposed 
by Metz and coworkers. Indeed the HC1CN complex can 
be formed at energies well below the abstraction reaction, 
however further reaction to form HC1 + CN would require 
surmounting barriers in excess of that for direct abstrac- 
tion. In order to explain the experimental results, which 
clearly show vibrational excitation of the CN, Harding 
proposes two possible mechanisms. First, the CN stretch- 
ing frequency is found in the calculations to be largely 
unaffected along the reaction coordinate. In contrast, the 
reagent CH stretch and product HC1 stretches, which are 
initially above the CN frequency, go through a minimum 
at the transition state and are below the CN frequency. 
Because these normal modes cross the CN stretching fre- 
quency, energy could flow between these two modes along 
the course of the reaction. Secondly, Harding's calcula- 
tions predict the existence of long-range entrance chan- 
nel complexes on both sides of the abstraction transition 
state (Cl-HCN, HCN-C1, and CN-HC1) which could also 
influence the final product state distributions. 

Most recently Troya et al. have performed a quasi- 
classical trajectory study of the CI + HCN — » HC1 + 
CN reaction |3lj . Their findings are in agreement with 
Harding, in that the intermediate HC1CN product only 
constitutes a minor channel (<5%) at the translational 
energies relevant to the previous experiments. The focus 
of the current study is to explore the long-range entrance 
channel complexes of X + HCN. 

EXPERIMENTAL 

A detailed description of the experimental apparatus 
has been given previously ^j|> therefore only the most 
salient features will be described here. Helium droplets 
are formed by expanding ultra-high-purity helium gas 
through a 5 //m nozzle which is cooled to 18-22 K. The 
helium stagnation pressure was maintained at 50 bar, re- 
sulting in the formation of droplets with a mean size of 
2500 - 6000 atoms j^j. Halogen radicals were produced 
by pyrolysis 0, m and doped into the droplets using an 
effusive source while the HCN is added to the droplets 
downstream using a simple scattering cell. The docu- 
mented growth of non-equilibrium structures in helium 
droplets is suggestive of a successive capture mechanism 
where dopants are pre-cooled before complexation occurs 
M- ESI • The relative timescales for species finding each 
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other in a droplet, to form a complex, and the rate of 
cooling are not precisely known. The cooling rate how- 
ever must at least be competitive with coagulation in or- 
der to prevent annealing of previously formed metastable 
structures in the pick-up process. In much the same way, 
the energy gained due to the mutual interaction of two 
dopants must also be quickly removed and complexation 
will occur at low interaction energies comparable to the 
droplet temperature of 0.4 K. Therefore, complexes are 
typically trapped in a local minimum reflecting the ap- 
proach geometry of the dopants; the low temperature 
does prevent rearrangement to the global minimum 

The infrared light from an F-center Laser (Burleigh 
FCL-20) or a PPLN-OPO (Linos Photonics OS-4000) 
interacts with the droplet beam using a linear multi- 
pass cell, designed to increase the effective interaction 
length. Excitation of the dopant inside the droplet leads 
to evaporation of approximately 600 helium atoms from 
the droplet 0], which reduces the on-axis beam flux 
reaching a liquid helium cooled bolometer. By mechani- 
cally chopping the laser and using phase-sensitive detec- 
tion, the resulting beam depletion can be recorded as a 
function of the laser frequency to obtain the absorption 
spectrum. A set of Stark electrodes was aligned orthog- 
onal to the laser interaction so that Stark and pendular 
spectroscopy could be performed. The DC electric field 
was oriented parallel to the laser polarization yielding 
AM=0 selection rules. Details on the calibration of the 
electric field are given elsewhere |45| . In order to optimize 
conditions for the pickup of a single halogen atom, the 
temperature of the pyrolysis source is adjusted, monitor- 
ing the percent dissociation of the precursor by probing 
both the X-HCN and X2-HCN binary complexes. 

ELECTRONIC STRUCTURE CALCULATIONS 

In order to characterize the precursor complexes we 
have calculated a two-dimensional angular potential (at 
the RMP2/aug-cc-pVDZ level) for HCN-C1 2 correspond- 
ing to the two rotors being fixed within the same plane. 
The two angles, 8 and <j>, were stepped in increments of 
10°, and the intermolccular distance R was relaxed to 
find its minimum in energy. The remaining geometric 
parameters were held fixed at the values found for sepa- 
rate HCN and Ck> optimizations, namely tch = 1-0779 
A, tq n~ 1.1828 A, and t C ici = 2.3275 A. The disso- 
ciation energy of the c omp lex was corrected for basis 
set superposition error |46J, and a spline interpolation 
was used to generate the final surface shown in Figure 
n The surface predicts a deep minimum for a linear ni- 
trogen bound geometry (6, <p) — (180, 180), and a very 
flat minimum around (90, 20) corresponding to a near 
T-shaped hydrogen bound complex. Previous gas-phase 
studies [ItI liil |49| have confirmed the linear geometry 
of the HCN-Cla, HCN-BrCl, and HCN-FC1 complexes 
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FIG. 1: A two-dimensional angular potential for HCN-CI2 
showing two possible local minimum structures. At each point 
only the intermolecular distance was optimized to find its min- 
imum in energy. Counterpoise correction was applied and the 
final surface generated using a bi-cubic spline interpolation 
algorithm. Contour lines are drawn in 100 cm -1 intervals. In 
helium droplets both isomers are observed for I2, Br2, and 
CI2 complexes with HCN, where as in the gas-phase only 
the global minimum (linear nitrogen bound) isomer has been 
identified |4?|. Interestingly only the global minimum HCN- 
F2 isomer was observed in helium. 

however there has been no previous evidence (nor predic- 
tion) of a hydrogen bound isomer. Fully relaxed geom- 
etry optimizations and harmonic vibrational frequency 
calculations confirm that the local minimum equilibrium 
structure is slightly bent away from an exactly T-shapcd 
geometry, however vibrational averaging is likely to be 
important due to the weak anisotropy. Exploratory cal- 
culations in which the restriction of planarity was relaxed 
did not result in any new minima. 

As pointed out by Harding 32], the reactivity toward 
HCN by halogen atoms is such that two reaction products 
have been postulated, with the halogen covalently bind- 
ing to either the carbon or nitrogen atom of HCN. Due to 
limitations of the basis set used and computational ex- 
pense, Harding admits that the calculations performed 
to date are more qualitative in nature. In this section we 
extend these calculations to higher-levels of theory and 
calculate each of the F, CI, and Br reactions so that they 
may be directly compared, but we only focus on the low- 
est energy channels including HXCN formation. To ex- 
plore the relative energetics of these complexes, we have 
employed the G2 method [5(J, and the resulting station- 
ary points, including zero-point energy, on the potential 
energy surface are shown graphically in Figure [21 For 
comparison, a few stationary points have also been cal- 
culated at the UMP2/aug-cc-pVTZ and UCCSD(T)/6- 
311++G(d,p) levels in addition to the G2 results, and 
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FIG. 2: Relevant stationary points of the X + HCN reactions calculated using the composite G2 method. Two reacted 
complexes in addition to two weakly bound entrance channel species are predicted. The HCNX products of chlorine and 
bromine atoms binding on the nitrogen end of HCN could not be stabilized in our calculations, however HCNC1 was found in 
reference [32^1 . The experimental heats of formation for the X + HCN — > HX + CN and X + HCN — > H + XCN reactions are 
taken from reference |51f . The photon energy used to study the pre-reactive X-HCN complexes is greater than the barriers to 
form HXCN, suggesting that it may be possible to photo-initiate this reaction. 



are summarized in Table [I] The calculated enthalpies of 
reaction are found to be in qualitative agreement with 
experiment at the G2 and UCCSD(T)/6-311++G(d,p) 
levels of theory, however the agreement is much worse at 
the UMP2 level. 

From Figure Hit is clear that the X + HCN -> HCNX 
reactions are all quite endothermic, and therefore we ex- 
pect that this reaction channel will not play at role in the 
dynamics at 0.37 K. More relevant to the experimental 
conditions is the reactivity of the halogen atoms with the 
carbon atom of HCN. For the F + HCN reaction, a very 
stable HFCN intermediate product is found in the calcu- 
lations. At the G2 level no barrier is predicted for this re- 
action, however at the UCCSD(T)/6-311++G(d,p) level 
the reaction barrier is 3.05 kcal mol -1 . Four of the six 
vibrational modes of this product were first observed by 
Andrews et al. [3^| in an argon matrix by codepositing F2 
and HCN, and then subjecting the matrix to broadband 
UV photolysis. Although rotational resolution was not 
achieved, the observation of the CN stretching vibration 
at 1672 cm -1 indicates that the CN bond has double- 
bond and not triple-bond character. Of interest to the 
current study, the CH stretching vibration was observed 
at 3016 cm -1 . Later Misochko et al. confirmed Andrews 
observations, and assigned the two remainin g vi brational 
modes which were not observed previously [3J|. Exper- 
imental evidence for a pre-reactive X+HCN complex or 
HC1CN however is still lacking. 



Harding's calculations on Cl+HCN are given in Table 
[I] and in general, agree well with our current calculations. 
G2 (UCCSD(T)) calculations predict that the reaction of 
a chlorine atom with the carbon of HCN is exothermic by 
-7.07 (-3.41) kcal mol -1 , and that the barrier height is 
3.30 (6.97) kcal mol -1 . For bromine atoms, the reaction 
to form HBrCN is predicted to be endothermic at the G2 
level of theory. Since the CH stretching vibration of the 
entrance channel complexes has a higher frequency than 
the barrier to form HXCN, it may be possible to photo- 
initiate the reaction starting from the entrance channel. 

While the above calculations give us an overview of 
the reaction energetics for the halogen atom - HCN re- 
actions, we are particularly interested in the topology of 
the surface around the entrance channel. For instance, 
it appears that HCN+X might be most reactive in a T- 
shaped arrangement, so the long-range orientational ef- 
fects of the potential may be very important in deter- 
mining the reaction dynamics. To begin to answer these 
questions we have calculated three-dimensional (R, 9, 
tch) non-relativistic adiabatic potential energy surfaces 
for Br + HCN and CI + HCN at the RCCSD(T)/aug- 
cc-pVDZ+{332} level of theory. The transformation of 
these surfaces to include spin-orbit coupling, and the cal- 
culation of the resulting bound states is being performed 
in collaboration with the group of van der Avoird, and 
details on these results will be presented elsewhere [52$ . 
In order to gain insight into our experimental results, we 
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G2 UCCSD(T) UMP2 Experimental Harding [3J 
6-311++G(d,p) aug-cc-pVTZ (AH) RHF+l+2+QC/cc-pVDZ 

CI (F) + HCN 0.0 0.0 0.0 0.0 0.0 

T.S. 3.30 (-1.70) 6.97 (3.05) - a - 8.2 

HC1CN (HFCN) -7.07 (-35.83) -3.41 (-28.44) 4.91 (-27.18) - -2.2 

HCNC1 (HCNF) - Q (15.32) - Q (22.60) - a (26.86) - 27.9 

HC1 (HF) + CN 24.21 (-9.45) 25.37 (-6.54) 42.90 (6.18) 21.41 (-11.58) 25.9 

C1CN (FCN) + H 23.11 (2.71) 26.80 (9.17) 19.03 (-3.03) 23.45 (9.26) 28^ 

"Could not be stabilized in the geometry optimization. 

TABLE I: A summary of the calculated energetics of the CI (F) + HCN reactions including zero-point energy. T.S. corresponds 
to the transition state of the reaction with the carbon atom of HCN to form HXCN. The experimentally determined heats of 
reaction (AH) are also shown for reference [sij . One can see that while the G2 and UCCSD(T) methods predict the enthalpies 
quite well, the UMP2 method shows much larger deviations. A graphical representation of these energies is shown in Figure [5] 
for the G2 calculations. Previous calculations by Harding for CI + HCN are also given for reference. 



will focus our attention on the qualitative features of a 
2D slice (r CH = 1.0655 A) from our 3D non-relativistic 
adiabatic potentials, which are shown in Figure [3] for Br 
+ HCN. The incorporation of the dimension correspond- 
ing to the CH bond length of HCN (tch) will lend insight 
into the change of the 2D potentials upon vibrational ex- 
citation of HCN (vi). 

Calculations were performed using MOLPRO 53] , em- 
ploying Jacobi coordinates, R and 9, and restricting the 
electronic wavefunction to C s symmetry, which gives rise 
to the 1A', 2A', and 1A" surfaces, corresponding to the 
three relative orientations of the unpaired orbital of the 
halogen atom with respect to HCN. A set of uncontracted 
mid-bond functions with exponents sp: 0.9, 0.3, 0.1, and 
d: 0.6, 0.2 (denoted as {332}) have been added to the 
aug-cc-pVDZ basis set and were placed at the midpoint 
between the halogen atom and the HCN nuclear center 
of mass [Hil IHij . For the plots shown in Figure |31 R and 
9 have been incremented by 0.1 A and 10° respectively. 
The ich and icn bond lengths of HCN were fixed at 
1.0655 and 1.1532 A respectively. Basis set superposition 
error has been accounted for by performing counterpoise 
correction 46] , and a spline interpolation used to smooth 
the surface. As noted by Fishchuk et al. for Cl-HF 56], 
there is some choice in how counterpoise correction is ap- 
plied due to the two A' states of the free halogen atom. 
Our procedure is the same as that reported for Cl-HF, 
in that the lowest RCCSD(T) energy of the free halo- 
gen atom was subtracted from both A' dimer energies, 
which preserves the double degeneracy of the n state for 
both linear geometries. The nitrogen bound (HCN-X) 
geometry corresponds to 9 = 0°. 

The surfaces clearly show that both linear geometries 
have stable minima, with well depths of -810 and -401 
cm -1 for the HCN-Br and Br-HCN isomers respectively. 
The isomerization barrier connecting Br-HCN to HCN- 
Br is found to be only 90 cm" 1 (neglecting zero-point ef- 
fects). At 6 = 0, the 1A" and 2 A' surfaces correlate with 
the doubly degenerate 2 H state, where as for 9 = 180, it 



is the 1A' and 1A" surfaces. The ground electronic state 
of the HCN-Br isomer is 2 £ and for Br-HCN it is 2 n in 
agreement with that found by Harding |32l] . These sym- 
metries illustrate that the unpaired orbital of the halogen 
atom is along the axis of the molecule for 2 E, while it is 
perpendicular to the molecular axis for 2 n. These differ- 
ent electronic symmetries are determined by electrostatic 
interactions of this dipole-quadrupole system. In the lin- 
ear nitrogen bound complex the negative end of the HCN 
dipole will want to approach the positive p-hole created 
by the unpaired orbital, thus orienting the unpaired or- 
bital on the molecular axis. When the positive end of 
the HCN dipole approaches the Br, the positive p-holc 
will be repelled, and orient itself perpendicular to the 
molecular axis, resulting in a n ground state. 

Fully relaxed geometry optimizations, at the RMP2 
level, were also undertaken for each of these stable struc- 
tures to generate dipole moments and scaled harmonic 
frequencies to compare with our experiment, which are 
summarized in Table ^ Note that vibrational frequen- 
cies could not be calculated for the basis sets which in- 
clude mid-bond functions, so we have instead used the 
aug-cc-pVTZ basis set. Calculations for iodine con- 
taining complexes employed the aug-cc-pVTZ-PP j5?J 
pseudo-potential for the iodine atom. One can see 
that while the binding energy for Br-HCN is in very 
good agreement with the 2D surfaces (399.5 cm -1 at 
the RMP2/aug-cc-pVTZ level compared with the 2D 
RCCSD(T)/aug-cc-pVDZ+{332} result of 401 cm- 1 ), 
the differences found for HCN-Br are somewhat greater 
(910 vs. 810 cm- 1 ). 

Relativistic - adiabatic potential energy surfaces 

The effects of spin-orbit coupling on reshap ing the po- 
tential energy surface is well documented [5^, and in 
order to make detailed comparisons with our experiment 
we must incorporate it. To include spin-orbit coupling 
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FIG. 3: The three non-relativistic adiabatic potential energy surfaces for HCN + Br calculated at the RCCSD(T)/aug-cc- 
pVDZ+{332} level. Contour lines are drawn in 50 cm -1 intervals. The surfaces predict both linear structures to be minima. 



RMP2 HCN-F HCN-C1 HCN-Br HCN-I (ECP) 

aug-cc-pVTZ (F-HCN) (Cl-HCN) (Br-HCN) (I-HCN (ECP)) 



D e (cm" 1 ) 


405.2 


834.3 


909.8 


935.1 




(189.0) 


(359.9) 


(399.5) 


(429.4) 


^harmonic (cill ) 


a 


3464.21 


3463.11 


a 




a 


(3443.85) 


(3433.52) 


(3430.88) 


^scaled (cm _1 ) 




3309.11 


3308.06 








(3291.66) 


(3281.78) 


(3279.26) 


R(A) 


3.240 


3.416 


3.542 


3.726 




(4.015) 


(4.342) 


(4.475) 


(4.670) 


B (cm" 1 ) 


0.1310 


0.0885 


0.0636 


0.0525 




(0.0881) 


(0.0561) 


(0.0405) 


(0.0339) 


A 4 (D) 


3.24 


3.66 


3.77 


3.91 




(3.18) 


(3.40) 


(3.47) 


(3.58) 



"Tight geometry convergence was achieved but imaginary frequencies were calculated 

TABLE II: Computed properties of the HCN-X (X-HCN) isomers at the RMP2/aug-cc-pVTZ level. The calculations for iodine 
incorporate a small core relativistic pseudo-potential (aug -cc-pVTZ-PP) in which 28 core electrons were replaced with the ECP 
[57)1 . Binding energies have been corrected for BSSE |4tl. Harmonic frequencies have been scaled by a factor of 0.9552, which 
is obtained by comparing a HCN calculation at this level to the corresponding experimental helium droplet band origin. R, is 
the X-HCN center of mass distance. 
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for the linear geometries (ID slices) of our 2D ab ini- 
tio potentials we take note of the operator form of the 
spin-orbit Hamiltonian: 



H so = AL-S 

= A(L Z S Z 



L+S- + L~S A 



) 



The spin-orbit coupling constant, A, is -269.3, -587.3, - 
2457, and -5068 cm -1 for F, CI, Br, and I, respectively 
[59l| . Note that the negative sign for A illustrates that the 
2 P 3 / 2 state lies below the 2 Pi/2 state in the free atoms. 
We make the approximation that the spin-orbit coupling 
is independent of the geometry of the complex, and re- 
mains that of the isolated atom. That this approxima- 
tion is valid for the long-range part of our potentials is 
justified because the electron-nuclei interaction scales as 
(^■)H3j an d therefore the inner lobes of the orbitals 
are more heavily weighted, which are expected to be 
only slightly perturbed by the relatively weak van der 
Waals interaction. Recent calculations on 0( 3 P) + HC1 
by Rode et al. 60] have examined the validity of this ap- 
proximation by explicit calculation of the spin-orbit ma- 
trix elements using the Breit-Pauli operator [6ll ]. Their 
results indicate that for internuclear distances similar to 
the van der Waals separation the spin-orbit coupling is 
not substantially affected. 

The diagonal term of the spin-orbit Hamiltonian 
(L Z S Z ) gives the 2 Hi/2 and 2 n 3 / 2 potentials as 2 I1 =F ^, 
while the 2 E potential is simply 2 E 1 / 2 - This is illustrated 
by plugging in the eigenvalues of the L z and S z opera- 
tors, namely A = 1 and E = +1/2 for the 2 n 3 / 2 state 
and A = and S = |1/2| for the 2 E 1/2 state. The off- 
diagonal spin-orbit term (L ± S ,=F ) couples the 2 II 1 / q state 
with the 2 £1/2 state and has the following form 
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To derive the relativistic adiabatic potentials we set up 
a R dependent 2x2 matrix for the SI = 1/2 states, set 
the diagonal elements to the energies from the ab initio 
calculation, and the off-diagonal elements to 2 1 / 2 ( : j), and 
diagonalize the matrix numerically to get its eigenvalues 
and eigenvectors. For this simple two level system, the 
potentials (eigenvalues) are given by: 



'E 1/2 = arn-^ + a-a'X'E) 
2 II 1/2 = a 2 ( 2 E) + (l-a 2 )( 2 H-|) 

where a is an R dependent mixing coefficient. The re- 
maining 2 n 3 / 2 state is simply 2 n + ^. The diagonal- 
ized relativistic and non-relativistic adiabatic potentials 
are shown in Figure At long range, where the E 
and II states are nearly degenerate, the 2 'E 1 / 2 state is 
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FIG. 4: Relativistic and non-relativistic adiabatic potential 
energy curves for the HCN-Br complex derived from the ab 
initio results. See text for details. 



pushed lower in energy compared to 2 E while the 2 n 1 / 2 
is pushed to higher energies compared to 2 n — ^4/2 due 
to the coupling. In the limit of R — ► 00 , the 2 n 3 / 2 
and 2 E]y 2 states converge and the separation between 
2 II 1 /2 and 2 n 3 / 2 becomes (3/2) A, or that observed in 
the free atom. The bound states of these potentials were 
calculated numerically using the program Level (6^ and 
the binding energies, rotational constants, and the van 
der Waals stretching frequencies are summarized in Ta- 
ble HTH To estimate the effects of CH stretching of HCN 
on the potentials, computations were also performed in 
which the HCN geometry was constrained to have a CH 
bond length of tch — 1-09 A, which is close to the vi- 
brationally averaged bond length for the v = 1 state. In 
general, the excitation of HCN results in a slight increase 
in the binding energy of the complex which we attribute 
to an increased dipole - quadrupole interaction. 

As noted above, spin-orbit coupling can strongly influ- 
ence the potentials, and indeed the dissociation energy 
of the HCN-Br ( 2 £i/2) complex is much smaller than in 
the non-relativistic calculation (383 vs 773 cm -1 ). The 
spin-orbit interaction is even strong enough to create a 
bound state in the upper 2 Ui/ 2 well, despite the fact 
that the 2 n state is purely repulsive. But since this state 
correlates with the excited spin-orbit component of the 
atom, we would not expect this isomer to be produced 
in helium. For the hydrogen bound isomer, the ground 
state is predicted to be 2 n 3 / 2 - Since this potential is not 
distorted from the non-relativistic 2 n, we expect that 
the nonrelativistic ab initio calculations for this isomer 
should be sufficient to compare with experiment. This 
may not be the case for the HCN-X complexes, and thus 
the nonrelativistic ab initio results should be interpreted 
with caution. After inclusion of the spin-orbit coupling, 
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Constant (cm -1 ) Br-HCN 
HCN vi = = 1) 2 n 2 n 1/2 2 E 1/2 

D 376.43 (390.04) 

B 0.03893 (0.03893) 

vvdws 47.37 (48.36) 



HCN-Br 



2 S 2 Si/ 2 2 IIi/2 

D 773.23 (775.89) 383.33 (384.94) 111.24 

B 0.06233 (0.06235) 0.05379 (0.05382) 0.04051 

u vdws 69.71 (69.84) 47.88 (47.97) 26.91 

a No bound states found 

TABLE III: A summary of the computed properties for the HCN-Br and Br-HCN isomers from bound state calculations on 
one- dimensional relativistic adiabatic potentials. Note that the computed bound states for the 2 n 3 / 2 state are the same as 
those for the non-relativistic 2 n state. 



we find that both isomers for bromine are now nearly 
isoenergetic. The smaller spin-orbit coupling constant 
for chlorine atoms allows the HCN-C1 complex to remain 
the global minimum. We look forward to the results of 
including spin-orbit coupling into our 2D and 3D sur- 
faces as it will be particularly interesting to investigate 
how the isomerization barrier is effected. Table HVI col- 
lects the molecular parameters derived from the bound 
states for the HCN-C1, HCN-Br, and HCN-I complexes. 

THE HCN-X 2 COMPLEXES 

Given that previous experimental 0, 0] studies have 
shown that halogen molecules bind to the nitrogen end of 
HCN, the search for the HCN-X2 complexes was straight- 
forward given the weak perturbation on the CH stretch- 
ing frequency. Pendular survey scans for Br2, CI2, and 
F2 + HCN are shown in FigureEl and the new peaks are 
labeled according to their gas dependence. The smooth 
variation in frequency shifts is telling of the magnitude 
of the interactions. Figure also shows the resulting 
zero-field spectra for the HCN-Br 2 , HCN-C1 2 , and HCN- 
F 2 binary complexes. Molecular iodine complexes with 
HCN were also observed but are not shown in the figure. 
Due to the fact that our calculations predict that each 
complex is linear, in accordance with the previous mi- 
crowave results |4g, |49j , the Q branches in the bromine 
and chlorine spectra are clearly unexpected. At this time 
we leave their explanation to a future study. We only pos- 
tulate here that an impurity complex could be overlap- 
ping the band, or that the helium interaction with these 
species is quite strong, due to the increased polarizability 
compared to Fa, and thus more strongly influences the 
spectra A similar effect in HCN trimer, another 

linear complex with a Q branch, has been treated theo- 
retically, and the Q branch was found to result from the 
thermal excitation of the droplet in the first few solva- 



tion shells, which rotates more rigidly about the a-axis 
[fiij . Possibly these molecular halogen complexes exhibit 
similar behavior. 

Longer pendular scans were also carried out to search 
the hydrogen bonded frequency region, given the predic- 
tion of a second, nearly T-shaped, minimum. Indeed, 
as shown in Figure for the case of bromine, a strong 
peak is observed at 3286.14 cm -1 which we assign to 
such a complex. Our assignment is based on pick-up 
cell pressure dependence measurements (optimized 1:1 
with HCN-Br2), scaled harmonic frequency calculations, 
rotational band contour analysis, and pyrolysis source 
temperature dependence. Hydrogen bonded complexes 
between molecular iodine and chlorine with HCN were 
also observed whereas the corresponding isomer of fluo- 
rine was not. The survey scans for all cases are shown 
in Figure as downward going peaks. The fact that a 
hydrogen bound F 2 -HCN isomer is not formed might be 
attributed to a smaller well depth, due to the smaller po- 
larizability of fluorine compared with chlorine, bromine, 
and iodine, and larger zero-point energy effects allowing 
it to convert back to the global minimum linear nitrogen 
bound geometry, which is observed. The infrared spec- 
tra for these complexes taken under field-free conditions 
will be presented elsewhere. The dotted lines in Figure 
illustrate the almost perfect linear scaling of the band 
origins for both the free stretch and the hydrogen bonded 
complexes. 

THE Br-HCN AND I-HCN COMPLEXES 
Br-HCN 

Given the identification of the precursor complexes, 
we are now in a position to heat the pyrolysis source to 
produce the corresponding atoms. As in our previous 
study of the X-HF systems, at sufficiently high pyrolysis 



Constant (cm 1 ) 


HCN- 35 C1 


HCN- 79 Br 




HCN-F 


v=0 (v=l) 


^ ^1/2 


2v 2v 

^ -^1/2 


2 E 


Sl/2 



Do 687.65 (689.98) 476.78 (478.86) 773.23 (775.89) 383.33 (384.94) 781.08 375.61 (377.24) 

Bo 0.08784 (0.08786) 0.08478 (0.08482) 0.06233 (0.06235) 0.05379 (0.05382) 0.05038 0.04245 (0.04247) 

Vvdw* 75.85 (75.98) 64.23 (63.53) 69.71 (69.84) 47.88 (47.97) 62.37 43.44 (43.52) 

A - -587.3 - -2457 - -5068 



a This calculation was performed with the aug-cc-pVDZPP+{332} basis set [£JJ. 



TABLE IV: A summary of the molecular parameters derived from the bound states of one-dimensional potential energy surfaces (see text). A is the atomic spin-orbit 
coupling constant used in the calculations and v v dms is the van der Waals stretching frequency, defined as the energy separation between the two lowest calculated 
bound states. For the "v = 1" calculations, the CH bond length was fixed at 1.09 A (compared to 1.0655 A for v = 0) in order to estimate the effects of vibrational 
excitation of the HCN. Note that the corresponding calculations for the X-HCN isomers are not presented since the ground electronic state of these complexes is 2 n 3 / 2 , 
which does not experience a first order coupling with the 2 E state, and so the molecular parameters are the same as those found in a standard ab initio calculation 
presented in Table IH1 
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FIG. 5: Pendular survey scans showing the formation of HCN-X2 molecular complexes (left panel) and their corresponding 
field-free spectra (right panel). Significant perturbations are observed in the HCN-Br2 and HCN-CI2 spectra (leading to the 
otherwise forbidden Q branch) which we ascribe to a helium droplet interaction. 
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FIG. 6: Pendular survey scans revealing the pyrolysis source 
temperature dependence (upward peaks are with a hot source, 
downward peaks are with a cold source) on the correspond- 
ing HCN + halogen experiment. The scans show peaks which 
we identify as molecular and atomic halogens complexed with 
HCN. The dotted lines are drawn to guide the eye, exempli- 
fying the linear scaling of the molecular halogen complex fre- 
quency shifts. The peaks labeled with an asterisk are known 
impurities. Vertical arrows point to the peaks which corre- 
spond to free CH stretch HCN-X complexes. 



source temperatures the signal levels associated with the 
molecular halogen complexes decrease, due to dissocia- 
tion into atoms, and a complementary set of peaks grow 
in. The pendular survey scans performed with a hot py- 
rolysis source are shown in Figure HO as upward going 
peaks. Each scan was taken at the appropriate pyrolysis 
source temperature for the dissociation of that particu- 
lar precursor. Concentrating on the bromine spectra, two 
new peaks are observed with the hot source, namely at 
3309.55 and 3284.61 cm" 1 . Based on the frequency shifts 
of these two bands it is already apparent that both linear 
isomers of HCN - Br may be formed, with the hydro- 
gen bonded isomer giving a much larger redshift on the 
CH stretching frequency than the corresponding nitrogen 
bonded isomer. We choose to first focus on the peak we 
preliminarily assign to a hydrogen bonded complex since 
a direct comparison can be made with the corresponding 
X-HF complexes studied previously [l2|. 

A field-free spectrum of the pendular peak centered at 
3284.61 cm -1 is shown in Figure[7fA). Again this band is 
only observed under conditions appropriate for bromine 
atom pickup. In agreement with that found for the Br- 
HF complex, P, Q, and R branches are observed, consis- 
tent with a linear complex in a n ground electronic state. 
For Br-HF a distinct 5B gap was observed between the 
Q branch transitions and those arising from the first P 
and R branch transitions which confirmed the ground 
state is 2 n 3 / 2 , consistent with the fact that the ground 
state of free bromine atoms is 2 P-$/2- The smaller rota- 
tional constant of Br-HCN however, makes these gaps in 
the spectra more obscure, but due to the predicted weak 
coupling of the spin-orbit interaction upon complexation, 



Constant 


Br-HCN 


I-HCN 


u (cm -1 ) 


3284.61(1) 


3277.79(1) 


B (cm" 1 ) 


0.0151(5) 


0.0120(5) 


D (cm" 1 ) 


1.5 x 10" 4 


1.0 x 10" 4 


[aA + (b+c)E] (cm" 1 ) 


0.04(1) 


0.04(1) 



TABLE V: A summary of the experimental parameters for 
the linear hydrogen bound X-HCN complexes obtained from 
a fit to the field-free spectra. 



we can say with confidence that this is also 2 n 3 / 2 band in 
analogy with Br-HF. A fit to the spectrum adopting the 
same Hamiltonian as used for Br-HF is shown in Figure 
EjB). A nuclear magnetic hyperfme perturbation of the 
form 

H' = aA(I ■ k) + b(I ■ S) + c(I ■ k){S ■ k) 

was also included in the Hamiltonian to reproduce the 
relative intensities of the P, Q, and R branches, in agree- 
ment with that observed for Br-HF, although here we 
do not resolve the individual hyperfme transitions. The 
(/ • k) term represents the coupling of the nuclear spin 
with the axial magnetic field of the molecule, which 
is proportional to A, (I ■ S) is the direct coupling of 
the nuclear and electronic spin angular momenta, and 
(/ • k)(S ■ k) is a second order interaction of the nuclear 
and electronic spins due to both of their projections on 
the molecular axis. Assuming Hund's case (a), the ex- 
pression for the interaction energy is: 



W = [aA+ (6 + c)S] 



9. 



J(J+1) 



I- J 



where 



I- J 



F(F + 1)- J(J + !)-/(/ +1) 



Note that the quantum number representing nuclear spin 
(I) equals 3/2 for bromine atoms, and F (the quantum 
number representing the total angular momentum in- 
cluding nuclear spin) is the coupling of J and I. The rela- 
tive transition intensities were calculated using spherical 
tensor matrix elements similar to those already reported 
for the Stark effect for X-HF [3 El, 

and a Boltzmann 

distribution for the rotational state population at 0.37 K. 
Individual transitions are convoluted with a 0.024 cm" 1 
Lorentzian lineshape function which was also varied to 
obtain the best fit. The spectroscopic constants (v : B, 
D, and [aA+ (b+c)S]) resulting from the fit are summa- 
rized in Table |V| 

The experimental (helium) rotational constant is a fac- 
tor of 2.7 smaller than that predicted from our ab initio 
calculations due to the fact that some of the helium fol- 
lows the rotational motion of the molecule, thus adding 
to the complexes moment of inertial |67j | . This magnitude 
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FIG. 7: (A) The rotationally resolved spectrum of the hydrogen bound Br-HCN isomer recorded under electric field-free 
conditions. The band shape is consistent with a 2 n 3 / 2 ground electronic state which is in agreement with theoretical predictions. 
The simulation (B) includes a nuclear magnetic hyperfine interaction (1=3/2 for bromine) due to the large magnetic moment 
of the bromine nucleus. A Lorentzian linewidth of 0.024 cm - was used for the simulation. 



of reduction is in good agreement with that observed for 
many other systems. The fact that the rotational con- 
stant of Br-HCN is reduced by a factor of 2.7 while the 
rotational constant of Br-HF is reduced by a factor of 
2.2 is also in agreement that the rotational constant re- 
duction factor is usually proportional to the magnitude 
of the B value in this range |6j|. For Cl-HF, the experi- 
mental vibrational frequencies have provided a very pre- 
cise benchmark for theoretical calculations [Hrl l6s(. and 
it will be interesting to compare similar calculations for 
X-HCN [j^. The scaled harmonic vibrational frequency 
(RMP2/aug-cc-pVTZ) for this isomer is 3281.78 cm" 1 , 
in good agreement with that observed (3284.61 cm -1 ). 

The vibrational frequencies of helium solvated rotors 
are typically very close to the gas phase values, owing to 
the weak interactions of the solvent, thus in general al- 
lowing a direct comparison with theory. A systematic in- 
crease in the vibrational red-shifts of complexes exhibit- 
ing linear hydrogen bonds has been observed however 
|67| . While a fully quantitative theory for this observa- 
tion is still lacking, such an effect can be qualitatively 
understood by the fact that the helium will act to reduce 
the amplitude of vibrational averaging, therefore making 
the hydrogen bond more linear, and thus stronger, lead- 
ing to a greater frequency shift. A second contribution 
to the observed red-shift is also expected to arise from a 
polarization effect due to the transition dipole moment 
of HCN, giving rise to dipole - induced dipole interac- 
tions with the solvent, which further lowers the energy 
of excited state with respect to the ground state. By 
fitting the enhancement of the frequency shift (helium - 
gas phase origin (Y)) versus the absolute magnitude of 
the shift observed in helium (complex - monomer origin 



(X)) for a database of 13 complexes observed both in he- 
lium and the gas phase, we have developed an empirical 
correction factor for the influence of the helium on the 
"gas-phase" band origin (in cm -1 ): 

Y = 1.822 + 0.03655Y (1) 

The standard deviation of Y is 1.2 cm -1 . See Figure 21 of 
reference for a plot of the experimental data points. 
Correcting for the helium induced shift one arrives at a 
predicted gas-phase origin of 3287.40 cm -1 . Although 
this makes the agreement with the scaled harmonic fre- 
quency somewhat worse, it is probably within the error 
of the ab initio calculations, and we await an estimate for 
the vibrational frequency from bound state calculations. 

I-HCN 

A similar study was carried out for the I-HCN com- 
plex, the infrared field-free spectrum being shown in Fig- 
ure A). The fit to the spectrum includes the hyperfine 
effects discussed above, however, 1 = 5/2 for iodine. The 
molecular constants derived from the fit are given in Ta- 
ble Here again the experimental rotational constant 
is consistent with the ab initio calculations, correspond- 
ing to a ratio of 2.83. The frequency shift from HCN 
monomer is measured to be -33.41 cm -1 (estimated gas- 
phase value using Equation^ -30.37 cm -1 ), in excellent 
agreement with the calculations employing the ECP (- 
31.94 cm -1 ). The nuclear magnetic hyperfine constant 
used in the fit is the same as for Br-HCN, namely 0.04 
cm -1 . This is reasonable given that the magnetic mo- 
ments for the atoms are quite similar |69ll70|. 
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FIG. 8: (A) A field-free spectrum of the I-HCN complex in 
helium droplets. The fit to the spectrum (B), which includes 
the effects of nuclear hyperfine (1=5/2 for iodine), was used 
to derive the spectroscopic constants given in Table A 
Lorentzian linewidth of 0.03 cm -1 was used for the simulation. 

From the survey scans of Figure^ peaks corresponding 
to the analogous Cl-HCN and F-HCN complexes are not 
observed. This is in contrast to the free stretch region 
where a new peak is observed for the HCN + CI case, 
which will be discussed below. 



THE HCN-Br AND HCN-I COMPLEXES 
HCN-Br 

We now turn our attention to the second set of new 
pendular peaks that grow in under the appropriate con- 
ditions for halogen atom pickup, namely those that ex- 
hibit a much smaller red-shift from HCN monomer. In 
the HCN + Br pendular survey scan (Figure[SJ), a peak is 
observed at 3309.55 cm -1 which is red-shifted from HCN 
monomer by -1.65 cm , in qualitative agreement with 
the calculations for a nitrogen bound HCN-Br isomer (- 
3.14 cm" 1 at the RMP2/aug-cc-pVTZ level). Because 
the ground electronic state of this complex is predicted 
to be 2 £i/2, spin-orbit coupling could be important in 
this system, and thus harmonic frequency calculations 
(which are performed on the 2 £ surface) are only semi- 
quantitative. From our ID PES's which include spin- 
orbit coupling, we find that the binding energy of the 
complex is increased by 1.61 cm -1 upon "vibrational ex- 
citation" of the HCN. Neglecting dynamical coupling of 
the different dimensions of the PES, this increase in the 
binding energy upon excitation necessitates a red-shift 
in the vibrational frequency of equal value, in excellent 
agreement with that observed, namely -1.65 cm -1 . Such 




u u u u u u 
Jf WimTl"hnjmffniCr? TTtifTml 



FIG. 9: Br-HCN and HCN-Br pyrolysis source pressure de- 
pendence measurements which aid in the determination that 
both observed bands are due to 1:1 complexes. 



excellent agreement could be fortuitous, and it will be in- 
teresting to compare this estimate with the results from 
higher dimensional calculations. An analogous compar- 
ison between the vch=0 and Vch = 1 bound states for 
the non-relativistic curves gives a frequency shift of -2.66 
cm , which is in much poorer agreement with experi- 
ment. 

The pick-up of dopants by helium droplets is a sta- 
tistical process, and the signal intensity as a function of 
dopant pressure, a Pick-Up Cell (PUC) pressure depen- 
dence curve, is given by a poisson distribution weighted 
by the droplet size distribution for a given set of source 
conditions j^. PUC curves are a valuable characteri- 
zation tool because higher order clusters can easily be 
distinguished by their optimum pressure. For our py- 
rolysis measurements, the pick-up zone is rather ill de- 
fined, however the optimum molecular halogen pressure 
upstream of the pyrolysis region could be used to aid in 
the determination of this new band to a second isomer of 
bromine atoms complexed with HCN. A pressure gauge 
was added to the pyrolysis source and the resulting pres- 
sure dependence curves for the two peaks at 3309.54 and 
3284.55 cm -1 are shown in Figure |5J Both bands are 
found to have exactly the same bromine pressure depen- 
dence, confirming that this new band at 3309.54 cm -1 is 
also due to a 1:1 complex. 

The field-free and Stark spectra for this band are 
shown in Figure HUT A) and (C). The band shape is consis- 
tent with that of a linear rotor in a 2 £ electronic state, 
as predicted by our earlier electrostatic arguments. In 
the field-free spectrum there appears to be fine structure, 
most notably as a splitting in the low rotational transi- 
tions. We tentatively ascribe this splitting to an inter- 
action of the (relativistic) 2 Xi/2 and 2 n 1 / 2 states, which 
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FIG. 10: (Left panel) (A) An infrared field-free spectrum for the CH stretching vibration of the HCN-Br complex. Fine 
structure is observed which we preliminarily assign to a parity splitting induced by spin-orbit coupling of the 2 £i/2 and 2 Hi/2 
states of the complex. A simulation of the band is also shown (B) which includes the effects of spin-orbit coupling. (Right 
panel) (C) The corresponding Stark spectrum for the HCN-Br complex recorded at an electric field strength of 2.25 kV cm -1 . 
The peak marked with an asterix is assumed to be related to an impurity complex. Lorentzian linewidths of 0.025 and 0.01 
cm -1 were used for fitting the field free and Stark spectra respectively. 



are coupled by the off-diagonal spin-orbit perturbation. 
The simulation shown in Figure ITUI was generated using 
an effective Hamiltonian which includes the off-diagonal 
spin-orbit interaction derived from one-dimensional cal- 
culations, however the splitting is not readily visible on 
this scale [jlj. Based on preliminary calculations of the 
bound states on the 3D spin-orbit corrected diabatic po- 
tential energy surfaces, it appears that the magnitude 
of the parity splittings for HCN-Br are smaller in a he- 
lium droplet than in the predicted gas-phase spectrum 
[5^|. Interestingly this behavior would be opposite to 
that observed recently for NO, where an increase in the 
parity splitting was observed jT^. The spectroscopic 
constants derived from the fit to the HCN-Br spectrum 
are collected in Table IVII The rotational constant ob- 
tained from the simulation is 0.019 cm' 1 , which can be 
compared with the zero-point corrected (Bo) values de- 
rived from the bound states on both the non-relativistic 
( 2 S) and relativistic ( 2 £x/2) ID potential energy sur- 
faces, namely 0.0623 and 0.0538 cm" 1 respectively. As 
shown in the ID potentials, and reflected in the B values, 
the spin-orbit coupling shifts the potential minimum to 
slightly longer internuclear distance. Since a fully quan- 
titative theory for the reduction of rotational constants 
due to the helium is still lacking, we cannot definitively 
say which of these calculated B values is in better agree- 
ment with experiment, since both reduction factors are 
within the range typically observed. By comparing the 
B reductions for the HCN-C1, HCN-Br, and HCN-I com- 
plexes however we do gain some insight into the accuracy 
of our relativistic potentials; see details below. 

A Stark spectrum for HCN-Br is shown in Figure 



HUT Ch recorded at an applied electric field strength of 
2.25 kV cm" 1 . The simulation includes the effects of the 
e/f symmetry splitting, however due to the mixing of 
J levels, its incorporation only slightly modifies the pre- 
dicted spectrum. The peak marked with an asterix is 
ascribed to an impurity complex and likely contributes 
slightly to the congestion in the field-free spectrum. The 
ground and excited state dipole moments obtained from 
the simulation are 3.78 D, in excellent agreement with 
that found at the RMP2/aug-cc-pVTZ level (3.77 D). 



HCN-I 

Based on the pendular survey scans shown in FigureEl 
we observe a pendular peak at 3309.37 cm -1 , which we 
tentatively assign a free stretch complex between HCN 
and an iodine atom. The field-free spectrum for this 
band is shown in Figure lllf A) . Again the band shape 
is consistent with a molecule with no net orbital angu- 
lar momentum (a S state), in agreement with the re- 
sults for the HCN-Br complex. The observed frequency 
shift from HCN monomer is -1.83 cm" 1 , which agrees 
nicely with the predicted increase in the binding energy 
upon vibrational excitation on the relativistic potentials, 
namely 1.63 cm" 1 . The rotational constant derived from 
the simulation (B) is B = 0.016 cm" 1 which is consistent 
with the results for HCN-Br. The corresponding value 
of the rotational constant derived from the lowest bound 
state of the 2 £i/2 potential is Bo = 0.04247 cm" 1 . Figure 
lllf C) also shows a Stark spectrum for HCN-I. The lines 
in the Stark spectrum appear to be much narrower than 



15 




FIG. 11: The field- free (A) and Stark (C) infrared spectra of the CH stretching vibration of the HCN-I complex observed in 
helium droplets. The simulation of the field free spectrum (B) was used to determine the molecular constants and includes 
an off-diagonal spin-orbit perturbation, which however only slightly modifies the band shape at this resolution. Lorentzian 
linewidths of 0.025 and 0.01 cm" 1 were used for fitting the field free and Stark spectra respectively. 



in the corresponding field-free scan, which we attribute to 
considerable unresolved fine structure in field-free scan, 
which is decoupled by the field. Such fine structure is 
most likely caused by the spin-orbit coupling in this sys- 
tem, in analogy to that observed for HCN-Br. The elec- 
tric field strength used to record the Stark spectrum (C) 
was not accurately calibrated and we can only say that 
the dipole moment is consistent with the ab initio value 
of 3.91 D. 



THE HCN-C1 COMPLEX 

When flowing CI2 through the hot pyrolysis source, 
no bands were observed that could be assigned to a hy- 
drogen bound Cl-HCN complex. In the free stretch re- 
gion however, a new peak is observed at 3309.33 cm -1 in 
good agreement with the previously observed bands for 
the HCN-Br (3309.55 cm" 1 ) and HCN-I (3309.37 cm- 1 ) 
complexes. The field-free and Stark spectra for HCN-C1 
are shown in Figure lT2[ A) and (C) respectively. This new 
peak is only slightly shifted from the HCN-CI2 band, and 
thus the pyrolysis source was run at the highest possible 
temperatures to reduce the CI2 related signals. Unfor- 
tunately these conditions also lead to an increase in the 
amount of HC1 produced in the source, presumably from 
the reaction of chlorine atoms with the quartz pyrolysis 
tube, and the HCN-HC1 molecular complex is also ob- 
served. The assignment of this band to HCN-HC1 was 
confirmed by a separate experiment in which HC1 was 
intentionally introduced into the pick-up cell. Despite 
the overlap of two impurity complexes, we are able to re- 
produce all of the features of the spectrum, and extract 
the molecular parameters which are listed in Table IVT1 In 



the analysis of the spectrum we were not able to identify 
any fine structure that could be attributed to a spin-orbit 
interaction, which we tentatively assigned in the bromine 
and iodine complexes. This is consistent with the much 
smaller spin-orbit coupling constant for atomic chlorine 
compared to bromine or iodine. Indeed a simulation of 
the HCN-C1 spectra using the effective Hamiltonian de- 
veloped for the HCN-Br, shows no such splittings at the 
experimental resolution. 

The vibrational band origin of HCN-C1 is shifted 
from HCN monomer by -1.87 cm . The ab initio 
frequency calculations predict a shift of -2.09 cm -1 at 
the RMP2/aug-cc-pVTZ level (see Table HTJ) in good 
agreement with the experimental value. The shift pre- 
dicted from the change in the binding energy upon vi- 
brational excitation from ID potential energy surfaces is 
-2.08 cm -1 . This value is nearly identical to the non- 
relativistic results due to the relatively small spin-orbit 
coupling constant for chlorine atoms. Bound state calcu- 
lations predict that the isotope splitting between HCN- 
35 C1 and HCN- 37 C1 is approximately 0.001 cm" 1 , which 
is negligible at our resolution. In contrast, the two iso- 
topes of Cl-HF were found to have a much larger splitting 
(0.038 cm" 1 ), due to the fact that it is a hydrogen bound 
complex [l3,ISa|- 

The rotational constant deduced from the fit to the 
spectrum is 0.032 cm -1 . Because the spin-orbit coupling 
is small in this case, the predicted rotational constants 
from the bound states on the 2 £ and 2 £1/2 potentials 
are very similar, namely B = 0.08478 (0.08784) cm" 1 
for HCN- 35 C1 on the 2 £ 1/2 ( 2 S) potential. The reduc- 
tion of the rotational constant upon solvation is thus a 
factor of approximately 2.7. If we assume that the solva- 
tion behavior between the HCN-C1, HCN-Br, and HCN- 
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FIG. 12: (Left panel) (A) A rotationally resolved spectrum of the vch band of the HCN-C1 complex in helium nanodroplets. 
The spectrum also has features due to both the HCN-HC1 and HCN-Cb complexes which overlap in this frequency region. 
Trace (B) is a summation of the 3 bands. (Right panel) (C) An experimental Stark spectrum of the corresponding band taken 
in the presence of a 2.833 kV cm -1 electric field. (D) A simulation of the experimental spectrum based only on the molecular 
parameters for HCN-C1. The peaks marked with an asterix are due to the overlapping bands. Lorentzian linewidths of 0.028 
and 0.013 cm -1 were used for fitting the field free and Stark HCN-C1 spectra respectively. 



I complexes is the same, then this reduction factor also 
gives support to our earlier comparison of HCN-Br and 
HCN-I with the B values obtained from the spin-orbit 
corrected potentials. On these potentials the reduction 
factor was also determined to be 2.75, however if we had 
compared the experimental value to the B value from the 
non-relativistic potentials, the B reduction would have 
been 3.2. 

Figure IT27 C) also shows a Stark spectrum for HCN-C1 
recorded at an applied field strength of 2.833 kV cm -1 . 
The dipole moment determined from the fit is 3.0(2) D, 
which is significantly smaller than the ab initio value of 
3.66 D. It is interesting to note that the ab initio cal- 
culations are in better agreement for bromine than for 
chlorine what could be due to vibrational averaging over 
the low frequency bending and stretching modes of the 
complex. This would be more important for CI due to its 
lighter mass, reducing the dipole moment from its equi- 
librium value. Indeed a similar trend was also observed 
for the X-HF complexes. 

DI-RADICAL Br-HCCCN-Br COMPLEXES 

Given the observation of the HCN-Br, Br-HCN, HCN- 
I, and I-HCN complexes, it is interesting to consider the 
possibility of picking up two halogen atoms, to form a 
di-radical complex such as Br-HCN-Br, which utilizes 
the HCN as a stabilizer, preventing the recombination 
to form HCN + Br2. Such complexes would represent 
the beginnings of nano-scale radical solids which could 
be used as high-energy density materials |73j . Gordon et 



Constant HCN-C1 HCN-Br HCN-I 



v (cm ) 


3309.33 


3309.55 


3309.37 


B (cm" 1 ) 


0.032 


0.019 


0.016 


D (cm" 1 ) 


5.0 x 10~ 5 


1.2 x 10~ 5 


1.0 x 10~ 5 


P (D) 


3.0 


3.78 





TABLE VI: The experimental spectroscopic constants ob- 
tained from a fit to the infrared spectra for each of the pre- 
reactive HCN-X complexes. Fine structure is observed in the 
HCN-Br and HCN-I spectra, which is attributed to spin-orbit 
interaction. A more detailed theoretical treatment for this ef- 
fect is given elsewhere |7l|. 



al. proposed a similar mechanism was responsible for the 
stabilization of nitrogen atoms in a N2 matrix (of up to 
10% atomic concentration) frozen in bulk liquid helium 
0, El El El El El While we have yet to perform 
extensive searches for HCN, we did explore this possi- 
bility for Br atoms and cyanoacetylene (HCCCN) , where 
better signals are observed. FigurefHflshows a set of pen- 
dular survey scans covering the important frequency re- 
gion. It is important to note that scans lT3T A) - (C) were 
recorded by picking up the HCCCN first, but in H3f D) 
the order has been reversed, and we instead pickup from 
the pyrolysis source first. In Figure ITS! A), only HCCCN 
is added to the droplets, while in H3r B). Br2 is flowing 
through a room temperature pyrolysis source. The result 
of heating the pyrolysis to the appropriate temperatures 
for bromine atom pick-up are shown in Figures lT3T C) and 
(D). In good agreement with the results of X + HCN, we 
find two peaks at 3326.1 and 3302.4 cm -1 (labeled by a 
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3296 3300 3304 3324 3325 3326 3327 



wavenumber (cm 1 ) 



FIG. 13: A series of pendular scans for Br + HCCCN where the bromine pressure is intentionally high to facilitate forming 
larger clusters. In (A) only HCCCN is picked up by the droplets and in (B) Br2 is flowing through the cold pyrolysis source. 
Scans (C) and (D) were recorded with the pyrolysis source at the appropriate temperature for bromine atom production, but 
for (D) the bromine is picked up first where as in (C) the cyanoacetylene is picked up first. The two peaks marked with a f 
are assigned to the linear Br-HCCCN and HCCCN-Br complexes, respectively. The two peaks marked with an * are found to 
optimize at higher HCCCN pressure and therefore correspond to more than one HCCCN. The peak at 3297.65 cm -1 optimizes 
at the same HCCCN pressure as the 1:1 complexes however is found to optimize at higher bromine pressure, suggesting that 
it is a complex containing two bromine radicals. The different behavior based on pick-up order is expected, since in (D) two 
bromine atoms will likely recombine to form Br2 in the absence of the HCCCN. 



f in the figure) which we assign to the Br-HCCCN and 
HCCCN-Br complexes, based on their frequency shifts 
and signal strengths. The peaks labeled with an * were 
found to optimize at higher HCCCN pressure, so they 
correspond to complexes containing more than one HC- 
CCN. The peak at 3297.65 cm -1 is a good candidate 
for being Br-HCCCN-Br because it only appears in the 
spectrum when the HCCCN is picked up first. This de- 
pendence on pick-up order is to be expected because two 
bromine atoms will likely recombine to form Br2 in the 
absence of the molecular spacer. To aid in this prelimi- 
nary assignment we performed bromine pressure depen- 



dence measurements which are shown in Figure IT-fl 

This new peak is found to optimize at higher bromine 
pressures than those peaks assigned to the 1:1 complexes, 
and agrees well with a simulation of a dimer using a Pois- 
son distribution for the pick-up statistics. For these pick- 
up "cell" pressure dependence measurements the main 
chamber ion-gauge was used as the x-axis, however since 
this really monitors recombined Bi~2 after collisions with 
the walls, instead of the local pressure of Br atoms in 
the pick-up zone, the pick-up curve for the dimer does 
not correspond to twice the optimum monomer pressure. 
Although we cannot rule out that a Br-HCCCN-Br2, or 
similar complex might have such a pick-up cell pressure 
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FIG. 14: Bromine pressure dependencies for the peaks as- 
signed to the Br-HCCCN, HCCCN-Br, and Br-HCCCN-Br 
complexes. The fact that the peak at 3297.65 cm" 1 optimizes 
at higher bromine pressure than the 1:1 complexes lends sup- 
porting evidence to our claim that this is a di-radical complex. 
The smooth curves are the calculated pick-up probabilities ob- 
tained by fitting a Poisson distribution to the 1:1 curve, and 
then re-simulating the curve for the dimer. 

dependence, a double-resonance population transfer ex- 
periment [80( in which this new peak is pumped and the 
HCCCN-Br2 complex is recovered would be definitive. 



DISCUSSION 

Extensive pendular survey scans have been performed 
in order to find evidence for the entrance channel com- 
plexes of fluorine atoms with HCN. From the ab initio 
calculations presented above, it is clear that the reaction 
to form the HFCN product is quite exothermic, and ex- 
hibits only a small barrier, if any, to the reaction. The 
calculated energies of the HFCN product and its barrier 
relative to F + HCN are -34.05 and -2.45 kcal mol -1 
(-28.44 and 3.05 kcal mol -1 ), respectively, at the G2 
(UCCSD(T)/6-311++G(d,p)) level of theory. Two en- 
trance channel complexes (F-HCN and HCN-F) were also 
found in the G2 calculations with energies -0.55 and -1.58 
kcal mol -1 respectively. The unintuitive result that the 
transition state is calculated (at the G2 level) to be lower 
in energy than the corresponding van der Waals minima 
is simply an artifact of the calculation, and this small 
energy difference is within the estimated accuracy of the 
method. The fact that we do not observe a fluorine atom 
pre-reactive complex could be an indication that indeed 
the insertion reaction takes place, even under our ex- 
perimental conditions at 0.37 K. This conclusion is also 
partially supported by the argon matrix work of Andrews 



[38j and Misochko .39] where the insertion product has 
been identified. No pre-reactive complexes were observed 
in these works. Guided by their observed frequency of 
the CH stretching vibration of 3018 cm -1 , we searched 
for HFCN in helium. No signals were observed, however 
this frequency range is at the edge of the tuning range of 
the lasers used in this work. The oscillator strength for 
this vibration is also approximately 50 times weaker than 
the HCN fundamental, based on harmonic frequency cal- 
culations, and thus the argon matrix experiment has the 
big advantage of concentrating the product using longer 
deposition times. We hope to revisit this experiment us- 
ing newer technology OPO lasers, with a larger tuning 
range and higher output power |8lLl82fl . 

Due to the statistical nature of the pick-up process it 
should be possible to "tag" the HFCN product with an- 
other weakly bound molecule, which could then act as a 
stronger infrared chromophore. Indeed we have prelim- 
inary evidence that such a situation also arises in reac- 
tion of aluminum atoms with HCN and HF Shown 
in Figure El is a two-dimensional angular potential for 
HCN + HFCN calculated at the RMP2/aug-cc-pVDZ 
level, which shows that multiple van der Waals isomers 
may indeed exist. Since cluster formation in helium oc- 
curs sequentially and obeys Poisson statistics, these clus- 
ters should be formed when picking up HCN after the 
fluorine atom, given that the droplet is large enough to 
stabilize the HFCN. Based on the calculated exothermic- 
ity, the formation of HFCN would boil off approximately 
2500 helium atoms (assuming 5 cm -1 per atom 0]). 

The last remaining question in this study is the exis- 
tence of a hydrogen bound Cl-HCN complex in analogy 
with that observed for bromine and iodine atoms. Due 
to the fact that the isomerization barriers between the 
hydrogen and nitrogen bound isomers are so small, one 
possible explanation is that this isomer is not formed 
due to a zero-point energy effect, and that the complex 
simply converts to the global minimum HCN-C1 isomer. 
To make any further analysis of this possibility we must 
await the spin-orbit corrected 2D potentials [52J, to get 
a more realistic view of how the isomerization barrier is 
effected by spin-orbit coupling. A second possibility is 
that the reactivity is enhanced when the chlorine atom 
approaches the hydrogen end of HCN and it reacts, in 
analogy with that proposed for F + HCN. 

SUMMARY 

We have presented here a combined experimental and 
theoretical investigation of the halogen atom - HCN en- 
trance channel complexes in order to compare with the 
X-HF complexes studied previously. The potential en- 
ergy surfaces for X + HCN are found to be considerably 
more complex than for X+HX, due to the possibility of 
hydrogen atom abstraction and addition reactions. The 
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FIG. 15: A two-dimensional interaction potential calculated at the RMP2/aug-cc-pVDZ level for the insertion product HFCN 
with a second HCN molecule. The second HCN molecule was assumed to lie in the plane of HFCN due to the fact that no 
non-planar geometries were observed in fully relaxed optimizations. The two angles were stepped in increments of 15° and 
only the intermolecular distance was allowed to optimize in the calculation. Counterpoise correction was applied and a spline 
interpolation used to generate the final surface. Contour lines are drawn in 100 cm" 1 intervals. By tagging the F + HCN 
reaction product with a second HCN molecule, we regain a strong high-frequency CH stretch which could be then probed 
experimentally. 



lone pair on the nitrogen of HCN is also found to signifi- 
cantly alter the shape of the long-range potential energy 
surfaces, and for X + HCN, we find that the global min- 
imum is a linear nitrogen bound geometry in contrast to 
the hydrogen bound global minimum for X-HF. A lin- 
ear hydrogen bound X-HCN complex is also predicted 
to be stable, however the isomerization barrier back to 
the global minimum is small. Ab initio calculations pre- 
dict that for iodine and bromine significant barriers to 
the chemical reactions exist, and thus the entrance chan- 
nel complexes should be quite stable when formed in a 
helium droplet. Both of the linear van der Waals iso- 
mers predicted by theory are observed experimentally for 
bromine and iodine atoms. Interestingly a HCN-C1 iso- 
mer was also observed whereas Cl-HCN was not. Since 
HCN-C1 is predicted to be the global minimum on the 
van der Waals potential, it could be that the isomeriza- 
tion barrier is simply too small to prevent Cl-HCN from 
rearranging to HCN-C1. A second possibility is that the 
reactivity is enhanced on the hydrogen end of HCN, and 
Cl-HCN goes on to react to form HC1CN. No entrance 
channel complexes were observed in the fluorine exper- 
iments, and we preliminarily interpret this as a result 



of reaction. The intermediate reaction product HFCN 
has been recently observed in an argon matrix, how- 
ever, also in that work no pre-reactive complexes were 
observed. Preliminary searches for HFCN in this work 
in helium nanodroplets did not reveal this species. The 
photon energy used to study the CH stretching vibra- 
tion of the entrance channel complexes is greater than 
the predicted barriers, suggesting that photo-initiation 
of the corresponding reactions could be possible. The 
unique asymptotic degeneracy of the multiple electronic 
states in the region of the entrance and exit channels for 
free radical - molecule complexes is found to give rise to 
fine structure in the observed spectra for HCN-Br and 
HCN-I. 
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